Landau level mixing by full spin-orbit interactions 
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We study a two-dimensional electron gas in a perpendicular magnetic field in the presence of both 
Rashba and Dresselhaus spin-orbit interactions. Using a Bogoliubov transformation we are able to 
write an approximate formula for the Landau levels, thanks to the simpler form of the resulting 
Hamiltonian. The exact numerical calculation of the energy levels, is also made simpler by our 
formulation. The approximate formula and the exact numerical results show excellent agreement 
for typical semiconductors, especially at high magnetic fields. We also show how effective Zeeman 
coupling is modified by spin-orbit interactions. 
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I. INTRODUCTION 

The manipulation of the spin of charge carriers in 
semiconductors, spintronics, has attracted increasing in- 
terest in recent years. In the paradigmatic Datta-Das 
spin transistor fi the spin of the electron passing through 
the device is controlled by the Rashba spin-orbit (SO) 
interaction^ which in turn can be varied by the applica- 
tion of gate voltages. The Rashba interaction stems from 
the structural inversion asymmetry (SIA) introduced by 
a heterojunction or by surface or external fields. In semi- 
conductors with narrower energy gap (InGaAs, AlGaAs), 
this effect is expected to be stronger. It has been shown 
experimentally that the Rashba spin-orbit interaction 
can be modified up to 50% by external gate voltages. 3,4 

In addition to the Rashba coupling there is also 
a material-intrinsic Dresselhaus spin-orbit interaction. 
This originates from the bulk inversion asymmetry (BIA) 
of the crystal, and can be relatively large in semiconduc- 
tors like InSb/InAlSb^ 

Both Rashba and Dresselhaus interactions contribute 
to the spin-dependent splitting of the band structure of 
the host material, leading to dramatic spin-dependent 
phenomena for electrons or holes in semiconductors: Ef- 
fects on the charge and magnetic transport, spin re- 
laxation, and spin-Hall conductivity have been recently 
studiediSiSi^iiSiii A drift-diffusive transistor, in contrast 
to the ballistic Datta-Das device, has been proposed to 
be more robust by Schliemann et alH for the case when 
the two SO interactions have the same strength. The 
experimental observation of the spin-galvanic effect and 
weak localization effects have increased the interest in un- 
derstanding the interplay between different SO terms— 
Photocurrent measurements have been used to obtain the 
ratio of Rashba and Dresselhaus coefficients^ 4 Beautiful 
optical measurements of SO effects have also been per- 
formed recently in strained semiconductors— 

It is well known that the determination of the eigen- 
values and eigenstates of the system is crucial for the 
calculation of a number of important physical properties 
of the system. In the case of a two-dimensional electron 



gas (2DEG) in a perpendicular magnetic field in the pres- 
ence of both SO interactions, few solvable cases have been 
analyzed in the literature (we devote the next section to 
review them). However, a comprehensive description of 
the more general case, in which both Rashba and Dres- 
selhaus interactions are present with arbitrary strength, 
would be interesting. 

In this work, we address the problem of SO coupling 
effects on the Landau level structure of a 2DEG in a 
strong perpendicular field. We find an excellent approxi- 
mate expression for the spinor Landau levels in the most 
general case of different Rashba and Dresselhaus inter- 
actions. We compare our approach to numerical results, 
made easier and recursively exact by our formulation of 
the problem. We further study the effective Zeeman g- 
factor of the system, and spin-orbit coupling is found 
to enhance or suppress the Zeeman splitting, depend- 
ing crucially on material parameters and gate voltages. 
This behavior may be found useful in the characteriza- 
tion of spin-filter— and spin-polarized currents in two- 
dimensional systems— 

In the next section we review the exactly solvable cases 
of SO coupling in a field, where either Rashba or Dressel- 
haus coupling is present. The third section is dedicated to 
the general case in which both terms are present. Using 
a Bogoliubov transformation we transform the Hamilto- 
nian of the 2DEG in a perpendicular magnetic field, and 
in the presence of both Rashba and Dresselhaus terms, 
to the 2DEG with only an effective Rashba interacting 
term with modified strength. The specific form of the 
interacting term allows for the derivation of the numeri- 
cal exact level structure of the model. Our approximate 
result is obtained and shown to be increasingly accurate 
for higher magnetic fields or weaker SO interaction. The 
last section contains typical results for different materials 
and discussion. 



II. SOLVABLE CASES 

The Hamiltonian of 2D electrons with effective mass 
m and Zeeman coupling g in a perpendicular magnetic 
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field Bz is 
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H = + v 
2m 2m 
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where fig is the Bohr magneton, a z is the Pauli matrix, 
and P = p + -A is the kinetic momentum. It is well 
known that this Hamiltonian can be written as a simple 
harmonic oscillator describing the Landau levels: 



H a = hhj c (a) a + — ) + fiiL> c £<T z , 



(2) 



in which a = (P y + iP x ) / ' \j2mhw c , u> c = eB/mc and 
6 = —g[iB'mc/2he. 

In this system, both Rashba and the linear Dresselhaus 
interactions have a simple form given by 



H R = ^-{P y cj x -P x a y ) = hw c a 
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Here we have introduced the dimensionless parameters 
« = «o v /^and/3 = /3o V / ^. 

Three exactly solvable cases are known in the 
literature^ 

a) When (3 — 0, the Rashba term couples the up-spin 
level (0„_i,O) to the down-spin state (0, (j) n ). The exact 
eigenstates and eigenvalues of H = Ha + Hr are spinors 
given by 



Vn 
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E r J l = hto c n T ^ \/l + ^ 2 nh 2 uj 2 /5 2 , 
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in which 8 — Tno c (l + 2£). Here n > 1 and the eigenstate 
n = exists only for ip l n with 0q = 0. The mixing angle 
n is given by tan 29 n = 2y/nahiu c /5 and varies from zero 
(when a = 0) to ir/4 (for weak magnetic field or strong 
SO interaction). Each Landau level 4> n ii splits into two 
levels tp l n and ip^+i with the splitting gap A = E l n —E^ +1 . 

b) a = 0. This time the Dresselhaus term couples 
(</>n,0) only to the (0, 4> n -i) leading to two new spinor 
states of H = Hq + Hd- 
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Here <5 = tuo c (— 1 + 2£), tan26*„ = 2^fnfihw c j '5 and n > 1. 
The eigenstate n = exists only for tp n with 9o = 0. 

c) The case of a = ±/3 is in general solvable in the 
absence of Zeeman coupling by going to a rotated spin 
coordinate a x ± a y } 2 For a — (3 the SO term is sim- 
ply H so = (a + 0o)(P x + P y ){a x - <r y ). In the diagonal 
space of a x — a y , the eigenstates are the harmonic oscilla- 
tor levels with a spin-dependent shift of the momentum 



Px/y — ¥ Px/y i ( a o + Po)m/h. For each Landau level in 
the absence of Zeeman coupling, the two spin-states are 
degenerate, but all states are shifted by a constant value 
AE = —Two c {a + /3) 2 /4. Note that there is no mixing 
between the two (rotated) spin states for different levels. 

In the presence of Zeeman coupling there is no ex- 
act solution of the energy spectrum. In a high magnetic 
field, where the Zeeman and spin-orbit splitting are small 
compared to the Landau level splitting, we can use per- 
turbation theory near the a — (3 point, and up to second 
order of perturbation we get (n > 0) 



El = hoj c [n + 1/2 + £ H ^ — 

V 7 ^ l + 2£ l-2£ / 

When 6 = we are left with the degenerate states with 
a constant shift of —fuv c f3 2 . 



III. GENERAL CASE 

In what follows we consider the general case when both 
spin-orbit terms are present. Without loss of generality 
we assume that \/3\ < \a\. Adding the Rashba interaction 
to the 2DEG Landau level results in levels as in JI}. In 
this case, keeping a constant while increasing (3, it is 
expected that the level splitting will decrease and the 
system will return back to the levels of ©, when j3 — > a. 

In the presence of both SO interaction terms, we use a 
Bogoliubov transformation, and introduce new operators 



and c' by 



c=(aa + il3a))/^Ja 2 - f3 2 



(7) 



for the case \(3\ < \a\. [For \(3\ > \a\ the definition of c 
and cr has to be interchanged.] This transformation is 

just the rescaling of the spatial coordinates x — ► xJ 



and y —> y\/ The total Hamiltonian transforms to 



H = H q (lu) + H so + V. 



(8) 



In this formula Hq(u>) is the Hamiltonian of a simple 
harmonic oscillator with modified energy u> = ui c (a 2 + 
(3 2 )/(a 2 — (3 2 ), which depends only on the ratio (3/a. 
The spin-orbit term H so is 



hw c \J a 2 
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which has the form of only an effective Rashba interaction 
term. Finally the spin diagonal interacting term V is 



V = huiX(c 2 



C f2 )A 



(10) 



in which the perturbation parameter A = 2a(3/(a 2 + ft 2 ) 
depends only on the ratio [3/a. We treat (|1(J|I as a pertur- 
bation term. The advantage of using transformation Q 
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lies in the simple form of i|10|) which makes it possible to 
achieve exact numerical solutions and most importantly 
a very accurate analytical ansatz for the spectrum. 

The exact eigenvalues and eigenstates of Ho(ui) + H so , 
which we use as the basis for the perturbation expansion, 
are simply given by those in Q , but with properly scaled 
frequency uj and SO interaction term 
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Here n > 1 and the eigenstate n — exists only for 
ip l n with 6q = 0. Note how the mixing angle and level 
splitting are renormalized by the Dresselhaus interaction. 
For (3 — > we recover the results Q . On the other hand 
near the point a = (3 both the angle and level splitting 
approach zero, as expected. 

Turning to the V matrix elements, we obtain 

A<^ = «|V|^)~*„ 4 , n ,± a , Pi/j =r,l. (12) 

Using the Brillouin-Wigner perturbation expansion, 19 
the basic formula of the wave-function is 



1*2*. 
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in which the prime requires rij ^ rij and pi ^ pj and the 
exact energy E^. has to be calculated from 



(14) 



The convenience of this perturbation scheme for V is that 
the task of finding the wave functions can be avoided^ 
Let us start from the second order correction to Ef}. : 



A (2) (££) = A 



,JHPi 12 



2 V - ^ \ u riiTij 



^ Ei 



Pi ,Pj ' 



(15) 



in which = E n \ — En]. For the fourth order cor- 

rection we first calculate the irreducible term 
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By irreducible we mean that the intermediate states are 
different from the original state ipf*. . There is also an- 
other contribution to the fourth order energy correction 
due to the reducible term 
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This reducible term is the product of two irreducible 
terms with complicated coefficients in general. The sim- 
ple form of the interaction term V allows us to derive 
(fTTfl just by replacing the bare energy E^. by the renor- 
malized energy E% + A^ 2 ) inside Ijl5|l and expanding up 
to A 4 terms: 

A( 2 '(£»+A( 2 ') = a( 2 )k;)+a«(^)+o(a 6 ). (18) 

It can be shown in general that all reducible terms of 
a given order can be derived from the irreducible terms 
of the lower orders if we use the renormalized energy 
in them. That means that we need to calculate only 
the irreducible terms. Once calculated, we renormalize 
the energy and iterate until we get the desired accuracy. 
Moreover, the calculation of the irreducible terms is sim- 
plified significantly because the transitions to intermedi- 
ate states are highly restricted, due to the peculiar form 
of (|12|1 . For a given level n the energy can be calculated 
with accuracy \( N ~ n ~)/ 2 7 where N is the total number of 
levels included (A < 0/a < 1). 

Knowing the limits for (3 = (@J and a = f3 © of 
the spectrum, we can also propose an approximate for- 
mula for the energy levels, increasingly accurate at high 
magnetic field or weak SO interactions. Both uj and A de- 
pend on f3 /a and only H so in (JSJ depends on the values 
of a and (3 as well. Rcscaling back @ , we notice that the 
effect of V is to transform Hq(u>) + V into the simple har- 
monic oscillator, with modified energy u>/ 'yl — A 2 = lo c , 
over the whole range of < A < 1. The main effect of 
V is to rescale the energy in (|HJ, but this consideration 
misses the constant shift of the energy levels near the 
point (3 = a. To compensate, we add an ad hoc energy 
shift which gives the correct level spectrum at a = /3 and 
at the same time does not affect the correct limit of the 
formula Q. We thus obtain 

1 
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+hw c f3 2 



l + 2£ 1-2^ 



El 



hw c n - - \/5 2 + An(a 2 - /3 2 )H 2 u> 2 
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in which 5 = hw c {\ + 2£), n > for E l n and n > 1 for E r n . 
As we will see in the examples, the agreement of Ijl9|) 
with the exact numerical results is excellent, especially 
at strong magnetic fields, when a and (3 are smaller. 



IV. RESULTS AND DISCUSSION 

In Fig. n the ground state and the first two excited 
states of a 2DEG in InSb are plotted. The magnetic 
field is chosen {B = 2T) such that the effect of spin-orbit 
splitting is visible and at the same time let the reader ap- 
preciate the accuracy of the approximate result (|19|l . In 
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FIG. 1: (Color online) Ground state and first two excited 
states for InSb as a function of ao/Qj 1 " 1 at po — (left 
panels) and as a function of 0o/ao a t oto = a™ ax (right 
panels). The fixed parameters are a™ ax = 250 meV A, 
g/go = —51, m/m e = 0.014, B = 2T. Energy is in units 
of hu c = 16.5meV. 



FIG. 2: (Color online) Ground state and first two excited 
states in InSb as a function of /3o//3cT a:c at ao = (left pan- 
els) and as a function of ao//?o at /3o = f3(J lo ' x (right panels). 
The fixed parameters are f3™ ax = 250meV A, g/g = —51, 
m/m e = 0.014, B — 2T. Energy is in units of hw c — 16.5meV 



the left part of this figure, first the parameter oq varies 
from zero to maximum value of a™ aa: = 250meV A while 
the Dresselhaus parameter is zero. Here the exact spec- 
trum |@J is plotted versus the Rashba interaction, akin 
to an applied gate voltage. Then keeping ag constant, 
on the right panels of the figure, we increase the value 
of (3q from zero to /?o = cxq. Here the circles show the 
exact numerical result, while the solid lines correspond 
to the approximate solution (|19|) . Note that as (3 can not 
affect the ground state energy directly, the correction to 
this level starts as a 2 /? 2 , while for all other levels, the 
correction is of order /3 2 and higher. 

In Fig. |2 we first switch on the Dresselhaus term (left 
panels) and then keeping (3 constant, we increase the 
value of the Rashba interaction (right panels) . This time 
it is the correction to the first excited level that is of 
fourth order a 2 (3 2 . 

Figure [3] shows how the SO interaction modifies the 
effective Zeeman coupling. Here this parameter is de- 
fined by g cS /g = 2{E l - EQ/(p, B Bg) where E l and E[ 
are the first excited state and the ground state respec- 
tively, i.e, the lowest spin-split level pairs. In all cases, 
the Rashba term increases g e g while the Dresselhaus cou- 
pling decreases it. The reason is simple as we note that 
the </>of state is coupled to <pii through Rashba inter- 
action and pushed further down, resulting in increasing 
spin splitting of each Landau level. Contrary to this the 
Dresselhaus interaction reduces the energy gap between 
the two spin states of each Landau level. Although not 
realized for the parameters chosen here, when the Zee- 
man splitting is small enough (as in low magnetic field), 
the Dresselhaus term can dominate the Zeeman coupling 
and interchange the ground state </>o| with </>ox. This is 
also found in the spectrum of quantum dotsS 
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FIG. 3: (Color online) The effective Zeeman coupling of the 
lowest spin-split levels in InSb, corresponding to Fig. Q B = 
2T. 



In general the effective Zeeman coupling depends on 
the strength of the spin-orbit coupling, bare Zeeman cou- 
pling and magnetic field. In the upper panel of Fig. 0] we 
have plotted this quantity for different semiconductors. 
The value of the magnetic fields is fixed at B = 2T for 
all cases. The bare value of the Zeeman splitting and the 
maximum value of the spin-orbit term is different and 
chosen as in the literature The larger the spin-orbit 
interaction, the larger is the change in the g c s. Notice 
that this modulation of g e g via the spin orbit effect is in 
principle controllable, via external gate voltage (to con- 
trol cto). This behavior may be useful in experiments 
with spinsiii In the lower panel of Fig. 0] we plot g e g 
as a function of the magnetic field when the value of the 
Rashba and Dresselhaus interactions is fixed at a™ ax and 
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/3q = ot™ ax /4, respectively. By increasing the magnetic 
field the effective Zeeman coupling decreases because the 
ratio of the spin-orbit splitting to the Landau level sepa- 
ration decreases. 




Once the energy is known the wave function can 
be calculated via (|f 3(1 which in turn can be used to 
calculate transport parameters like charge and spin 
conductivityiii*22i2i That work is in progress. 
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FIG. 4: (Color online) Upper panel shows the effective Zee- 
man coupling for different semiconductors as a function of 
only ao /a™ ax (left) and as a function of (i /ct™ ax but in the 
presence of Rashba interaction. The lower figure is the be- 
havior of f/ c ff as a function of magnetic field for fixed value of 
spin-orbit interactions. 
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